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ABSTRACT 



Oh- 



Context. Based on their spectral energy distribution, Herbig stars have been categorized into two observational groups, reflecting their 
overall disk structure: group I members have disks with a high degree of flaring as opposed to their group II counterparts. Literature 
^0 ■ results show that the structure of the disk is a strong function of the disk mass in pm-sized dust grains. 

Aims. We investigate the 5-35 yum Spitzer IRS spectra of a sample of 13 group I sources and 20 group II sources. We focus on the 
continuum emission to study the underlying disk geometry. 

Methods. We have determined the [30/13.5] and [13.5/7] continuum flux ratios. The 7-fim flux excess with respect to the stellar 
photosphere is measured, as a marker for the strength of the near-IR emission produced by the hot inner disk. We have compared our 
5— i ' data to the spectra produced by self-consistent passive-disk models, for which the same quantities were derived. 

Results. We confirm the results by Meijer et al. (2008) that the differences in continuum emission between group I and II sources can 
largely be explained by a difference in amount of small dust grains. However, we report a strong correlation between the [30/13.5] 
and [13.5/7] flux ratios for Meeus group II sources. Moreover, the [30/13.5] flux ratio decreases with increasing 7-fim excess for all 
targets in the sample. To explain these correlations with the models, we need to introduce an artificial scaling factor for the inner 
disk height. In roughly 50% of the Herbig Ae/Be stars in our sample, the inner disk must be inflated by a factor 2 to 3 beyond what 
hydrostatic calculations predict. 

Conclusions. The total disk mass in small dust grains determines the degree of flaring. We conclude, however, that for any given disk 
mass in small dust grains, the shadowing of the outer (tens of AU) disk is determined by the scale height of the inner disk (~1 AU). 
The inner disk partially obscures the outer disk, reducing the disk surface temperature. Here, for the first time, we prove these effects 
observationally. 

Key words, circumstellar matter; Stars: pre-main sequence; planetary systems: protoplanetary disks 
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• • 1 . Introduction )millimet er emission are s ignificantly smaller in group I than in 

. ^ ; group II (lAcke et al.ll2004l) . 

K,* Herbig Ae/Be stars are the intermediate-mass (2-8 M B ) ana- T ,. , , , ... , , , 

rN i c rr, rw, • . rp, a a u ■ »i In the early disk models, the inner disk was modeled as a 

s_] logues or T Taun stars. They are surrounded by circumstel- 

53 . i j- i u- u u v j 1 u lL c ■ i * vertical wall at the dust sublimation distance to the star. Later 

CZ lar disks, which are believed to be the sites or ongoing planet „ . it — - — — — i u»»a . . , 

< c .■ rp, . e i - i i i u- * refinements (e.g., Isella & Natta 2005) took into account the de- 
formation. The geometry of these disks has been the subject , . v 6 ." — : — 531 " 

c . .. rri TT1 /PtTwvTK u i c a \u pendence of grain sublimation temperature on gas density, re- 

of numerous studies. Meeus et al. (2001) have classified the r , . . b , . . , r , 

, j . , j. . ,n rr , , c if suiting in a curved inner rim. With the coming of age of near 

observed spectral energy distributions (SEDs) of a sample of T _ . b „ , . , , , , 

u , . . • . . t i j. i ■ IR interferometers, however, it became clear that our knowl- 

Herbig stars into two groups: group I members display a sig- 



..c , c ■ c a Vu .u • tt * edge of the inner disk remains incomplete. For a growing list 

^^EM^^ ^Z^Z 1 ! ofWs (Monmer et al. 2005; Kraus et al. 2008- Acke et al. 



. — " — , p p. r-rn u nr . A i , a a u 1 01 targets Monmer et at. 2UU5; Kraus et al. 2UU8; Acke et al 

parts. Dullemond & Domimk (2004a) have provided a physical „„„„ fc — ■ . , , II A^r.ol I — : T~\ — \ 

, . • e .,• j-rr tt u * 2008 ; Tanmrkulam et al. 2008 ), the inner disk appears less spa 

explanation for this difference: group II sources have an outer 1 1 1 



j. , , . , . . . j • ' j. . . ,. .. , tially resolved than expected from the models. This is corn- 
disk which is protected against direct stellar radiation by the \ , • , , , r 

cc , j- i tt -c ,u a- i c monly explained by the presence of an optically thick, possibly 

pufted-up inner disk. However, if the outer disk emerges from J \. , 1 r , 

f, • j- , , , ■ • , , a ■ , ., (, cr , gaseous, disk component located within the dust sublimation ra- 
the inner disk s shadow, i.e. has a large flaring angle, its SED re- IT TT, • f • , 

. c T ip. ii ao r>. — ^O/ lonn/iuk dius. However, near-IR interferometric data have been success- 

sembles that of a group I source. Dullemond & Domimk (2004b) _ „ „ . . , . n . . ,. 

j| A/r ■■ m Tonno K ■ a ' *■ +v. + +v. a » • • a- * u fully fitted with artificially increased inner disk heights as well 

and Mener et al. 2008) indicate that the dust grain size distribu- nT J . , , • , . „ . ^ 

.. ' T~ ■ ■ , • t . . . ftU (Verhoeff et al. 2009, submitted to A&A). 

tion plays a determining role in the structure of these circum- v ' 

stellar disks and their SED appearance. This idea is supported We have col lected and analyzed the Spitzer IRS 

by observational evidence: the grains responsible for the (sub- dHouck et al.l [2004) spectra of a sample of 54 Herbig Ae/Be 

stars (Bouwman, Juhasz, Acke et al. in prep). In this Letter, we 

focus on the shape of the underlying continuum and report the 

* Postdoctoral Fellow of the Fund for Scientific Research, Flanders. discovery of two tight correlations between the near-IR excess 
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HD244604 




15 20 25 
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Fig. 1. Spitzer IRS spectrum of the sample target HD244604. 
Indicated are the continuum wavelengths at 7, 13.5 and 30 fim 
where the spectral features of silicates and PAHs are weak or 
absent. The bin width is 1 fim. 



and mid-to-far-IR colors. We explain these correlati ons in the 
frame work of self-consistent passive-disk models (Mi n et ail 
2009). 



2. Anti-correlation between near-IR excess and 
far-IR color 



For all the target stars, we have fitted a r eddened K uruczl dl99ll) 
model to optical photometric measurements available in the lit- 
erature. From the IRAS 60 fim photometry, the excess mag- 
nitude is computed (M exc = 2.5 log(F tot /F*)). The Spitzer IRS 
spectra are used to determine the average fluxes at 7, 13.5 and 
30 fim in a bin with a width of 1 fim. As shown in Fig. Q] only 
a few wavelengths are suitably free of solid-state features. The 
chosen bin width is a compromise between maximal signal-to- 
noise and minimal influence of spectral features of amorphous 
and crystalline silicates, and PAHs. The error bars on the mea- 
sured fluxes include the errors in the individual spectral channels 
and the spread of the spectrum within the bin. In this paper we 
discuss the [13.5/7] and [30/13.5] flux ratios, and the 7-fim ex- 
cess which probes the near-IR excesfl The error on the latter 
includes an addition al 0.1 1 mag, as the ph otometric accuracy of 
Spitzer IRS is 10% dBouwman et al.ll2008l). 

In previous work dAcke & van den Anckerl l2006), we have 
shown that the IRAS 60 fim exces s can distinguish in a simple 
and reliable way between the two iMeeus et al.l d2001l) groups. 
Targets with a 60 fim excess exceeding 10. 1 magnitudes are clas- 
sified as group I, the others as group II. Seven sources with tenu- 
ous or optically thin disks, defined as those with a total infrared- 
to-stellar luminosity Lir/L* below 10%, are excluded as they 
cannot be compared to our optically thick disk models (Sect. [3]). 
Our final sample contains 13 group I and 20 group II sources for 
which Spitzer IRS spectra are available covering the full range 
from 5 to 35 fim. 

In Fig. [2]3, the [30/13.5] flux ratio is plotted as a function of 
the excess at 7 fim. For both group I and group II sources, a clear 
decline in [30/13.5] flux ratio is observed with increasing 1-fim 
excess. As indicated by the Pearson correlation coefficient, the 
trend is significant with a false-positive probability (FPP) less 



Table 1. Model parameter ranges, m is the mass in sub-fim- 
sized dust grains, p is the power which governs the decline of the 
surface density with increasing distance to the star (S oc R~ p ), i 
is the inclination of the disk with respect to the line-of-sight, <I> 
is the artificial scaling factor for the inner disk scale height, with 
<J> = 1 indicating hydrostatic-equilibrium models. See text for 
details. 



Parameter 



Values 



Approx. 70% of the 7-fim flux is produced in the innermost disk 



logra[M ] -3,-4,-5,-6,-7 

p 1,1.5,2 

i 10° - 90° 

O 1,2,3 



than 1%. We have performed a linear fit to the {1-fim excess, 
log[30/13.5]) data of group I and II separately. The fits, indi- 
cated by the black and red line in the figure, have similar slopes: 
-0.16 + 0.01 and -0.18 + 0.01 for group I and II respectively. 
Note, however, that the [30/13.5] flux ratios in group I are higher 
than those in group II for any given value of the 1-fim excess. 

Fig. |3? is a plot of the [30/13.5] versus [13.5/7] flux ratio. 
Group II sources display a pronounced correlation between both 
colors (FPP < 1%). The group I sources are scattered in this 
diagram. A few outliers are noted in the plot. Group I sources 
HD34282 (#2) and HD135344B (#10) are the bluest sample 
sources in [13.5/7] color, and the reddest in [30/13.5]. These 
targets have indications for the presence of a large (opacity) 
gap in their inner disk, which explains their anomalous colors. 
HD3741 1 (#20) and HD 152404 (AK Sco, #29) have the highest 
[30/13.5] ratio in group II and deviate from the trend. Also for 
these sources, a large disk gap may be present. The presence of 
a (stellar or sub-stellar) binary companion likely plays a role in 
opening up the gap. 



3. Comparison with disk models 

From the correlations described above, it appears that the hot 
inner disk (~1 AU) and the structure of the outer disk (tens of 
AU) closely relate to each other in the sample of Herbig stars. 
To understand the physics behind this result, we have computed 
a grid of disk m odels using the radiative transfer code MCMax 
(Min et al. 2009). The central star in these models is a main- 
sequence star of spectral type AO (T e ff = 10000K, R* = 2 Rq, 
M* = 2.5 M ), representative for the median star in the sam- 
ple. Altering the stellar parameters over the sample range has no 
influence on the analysis below. 

A circumstellar disk was modeled, assuming hydrostatic 
equilibrium (HSE) and thermal coupling between gas and dust. 
The gas-to-dust mass ratio was fixed to 100. The dust consists 
of astronomical silicates with a grain size of 0.12 fim. Large (> 
several fim) grains have not been included , as they do not in - 
fluence the SED shape shortward of 60 fim (Me iier et alj |2008). 
The total dust mass in the models (m) hence only refers to the 
total dust mass in these small grains and should be interpreted 
as such. We have computed models with different values for m. 
Next to the dust mass, the surface density power law (E oc R~ p ) 
and the inclination of the system (/) were altered. Also, we have 
artificially increased the inner disk's scale height by a factor O, 
to study the effect of a more puffed-up inner rim than HSE calcu- 
lations produce. The ranges of the model parameters are listed in 
Table [1] As for the sample targets, only disk models that display 
a 10 fim silicate emission feature in their spectrum are included 
in the analysis. In total, 1000 disk models were included in the 
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Fig. 2. Left panels: (a) The clear anti-correlation between the 1-fim excess and the Spitzer [30/13.5] flux ratio. Meeus group I and 
II sources are resp. black and red dots. Also shown are the best linear fits to the correlation for group I (black) and II (red line). 
Group I: 1. HD31293, 2. HD34282, 3. HD36112, 4. RR Tau, 5. HD38120, 6. HD250550, 7. HD97048, 8. HD100453, 9. SS73 
44, 10. HD135344B, 11. HD139614, 12. HD142527, 13. HD169142, Group II: 14. HD31648, 15. HD35187, 16. HD244604, 17. 
HD37258, 18. BF Ori, 19. HD37357, 20. HD37411, 21. HD37806, 22. HD72106S, 23. HD85567, 24. HD95881, 25. HD101412, 
26. HD104237, 27. HD142666, 28. HD144432, 29. HD152404, 30. VV Ser, 31. WW Vul, 32. HD190073, 33. HD203024. (c) and 
(e): Identical to panel (a) with all models included (dots). The stars refer to the observations, (c) The color of the dots refers to 
the total disk mass in small dust grains: log m [M ] = -3 to -7 from red to blue. The models shift from upper right to lower left 
with decreasing mass in small dust grains, (e) Colors refer to the inner disk scale height parameter <E>: yellow are the models in 
hydrostatic equilibrium, green and blue refer to O equal to 2 and 3. The models shift from upper left to lower right with increasing 
inner rim scale height. Right panels: (b) The [30/13.5] flux ratio as a function of the [13.5/7] ratio. A tight correlation between both 
flux ratios is noted for group II. (d) and (f): Identical to panel (b) with models included (dots). The same plotting convention as in 
the left panels is used. The models shift from upper left to lower right with decreasing mass in small dust grains (d) and from upper 
right to lower left with increasing scale height (f). 



study. Based on their 60 /im excess, roughly half would be clas- 
sified as Meeus group I sources, the other half as group II. 

Varying the inclination of the considered disk models only 
has a minor influence on the 7-yum excess (variations of 0.2 mag) 
and the IR flux ratios (25%). A spread in inclinations there- 
fore cannot reproduce the range, location and correlations of the 
observations in the diagrams. The locus of the models is pre- 
dominantly determined by the mass in small dust grains (see 



Fig. Wc>d). The higher the mass in small dust grains, the larger 
the outer disk's flaring angle becomes, and hence the redder the 
model appears in the [30/13.5] flux ratio. However, the dust mass 
cannot account for the correlation that is seen between the 7-yum 
excess and the [30/13.5] flux ratio in both groups, and the corre- 
lation between the [13.5/7] and [30/13.5] flux ratios in group II. 

The power of the surface density law has a moderate in- 
fluence on the location of the models in the diagrams: models 
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Fig. 3. The distance to the best linear fit to the group II data in 
Fig. [SJj (red line) is strongly correlated with the 60-pm excess. 
The black line represents the best fit to the group I data in the 
same figure. The stars are the group I (black) and group II data 
(red), the dots are the models. Colors reflect the mass in small 
dust grains, from m = 10" 6 (bottom left) to 10~ 3 M (top right). 



with a steep decline of the surface density with radius are bluer. 
However, the changes induced by this parameter are insufficient 
to explain the observed spread in disk excess and color. To cover 
the full observational range, especially in 7-pm excess, we need 
to artificially inflate the inner disk beyond its equilibrium scale 
height by a factor of a few (see Fig. [2k/). 

The correlation between the [30/13.5] color and the 7-pm 
excess displays scatter that greatly exceeds the observational er- 
rors. We take the best linear fit to the group II data (i.e. the red 
line in Fig. [2?) as the reference and compute the perpendicular 
distance of all observations to this line. This distance is strongly 
correlated with the magnitude of the 60 pm excess: sources be- 
low the average trend have lower excesses than targets above 
it (see Fig. [3). This correlation with far-IR excess is significant 
for the sample as a whole (FPP < 1%), but also when limited 
to group II sources alone (FPP < 5%). The disk models also dis- 
play the increase in 60 /mi excess with distance to the average fit. 
This confirms that the mass in small dust grains determines the 
locus of the targets and the disk models in the color-excess dia- 
gram, perpendicular to the observed correlation. Group I disks 
are more massive in small dust grains (m > 10- 5 M ), while 
the group II disks cluster around m = 10~ 6 - 10~ 5 M©. Along the 
color-excess and color-color correlations, the defining parameter 
is the scale height of the inner disk. Roughly half of the targets 
can be modeled with a disk in hydrostatic equilibrium. The disks 
that are blue in [13.5/7] and [30/13.5] colors, however, have large 
7-pm excesses. For these sources, an artificially increased scale 
height of the inner disk is required in our models. 



4. Discussion 

The degree of flaring is the increase in opening angle of the disk 
surface a(r) (= arctan(/;(r)/r), with h(r) the height of the surface 
above the midplane) with increasing distance to the star r. The 
disk surface is the surface that radiates in the infrared, where the 
opacity along the line-of-sight reaches unity. All the model disks 
flare, but the flaring is stronger in models with a high mass in 
small dust grains (h{r) oc r 130 vs. r 115 ). Moreover, the maximal 



opening angle of a massive disk is approximately 5 times that of 
a low-mass disk (a max = 25° vs. 5°). 

The observed range in 7-pm excess surpasses the range pre- 
dicted by hydrostatic-equilibrium models. An inflation of the 
inner disk height, beyond what HSE models produce, naturally 
leads to a higher near-IR excess, but also affects the outer disk. 
A disk with a puffed-up inner rim casts a large shadow over the 
outer disk. The surface layers of the outer disk are therefore cool 
and will produce hardly any infrared excess in the mid-to-far- 
IR. The combination of both effects reduces the [30/13.5] and 
[13.5/7] flux ratios with increasing 7-pm excess and explains the 
observed correlations. Note that the opening angle and degree of 
flaring of the disk is set by the mass in small dust grains, and 
only marginally depends on the height of the inner disk: a max 
decreases by 10% when <t) is increased from 1 to 3 for a given 
model. 

5. Conclusions 

Using Spitzer IRS spectra of a large sample of Herbig Ae/Be 
stars together with a large parameter set of radiative transfer 
models we come to the following conclusions. 

- There is a strong observational correlation between the near- 
infrared excess and the infrared spectral slope of Herbig 
Ae/Be stars. This indicates that the apparent geometry of the 
outer disk is determined by that of the inner disk. 

- Using radiative transfer modeling we show that the structure 
of the disk is governed by two parameters: 

1 . The mass in small dust grains which determines the de- 
gree of flaring. 

2. The inner disk scale height which determines the degree 
of shadowing. 

- Disk models in hydrostatic equilibrium fail to reproduce the 
observed spread in inner disk scale heights. 

The findings presented in this letter together with literature 
infrared interferometric results, show that our understanding of 
the innermost regions of Herbig stars is yet incomplete. It is clear 
that additional physical processes, not included in our models, 
play a role. The correlation between the 7 /mi excess and the 
apparent geometry of the outer disk shows that the mechanism 
causing the flux at short wavelengths influences the outer disk. 
This rules out optically thin hot gas emission as an important 
contributor to the 7 pm excess in Herbig stars. Interestingly, re- 
cent modeling effo rts including gas chemistry calculations (e.g. 
Woitke et al. 2009) show that the gas temperatures can be much 
higher than the dust temperatures in the inner disk. A similar ef- 
fect is obtained in disk model s when including the stell ar far-UV 
to X-ray radiation field (e.g. lGorti & Hollenb ach 2009). As a re- 
sult the inner disk height is strongly increased, while the outer 
disk remains unaltered. Further study is needed to verify whether 
this affects the disk continuum emission in conformance with the 
observations. 
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